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ABSTRACT 



We report on the abundances of 13 elements in the planetary host HD 13189, a 
massive giant star. Abundances are found to be sub-solar, with [Fe/H] = — 0.58± 
0.04; HD 13189 is one of the most metal-poor planetary hosts yet discovered. 
Abundance ratios relative to Fe show no peculiarities with respect to random field 
stars. A census of metallicities of the seven currently known planet-harboring 
giants results in a distribution that is more metal-poor than the well-known 
metal-rich distribution of main sequence (MS) planetary hosts. This finding 
is discussed in terms of accretion of H-depleted material, one of the possible 
mechanisms responsible for the high-metallicity distribution of MS stars with 
planets. We estimate the mass of the HD 13189 progenitor to be 3.5 Mq but 
cannot constrain this value to better than 2-6 Mq. A stellar mass of 3.5 Mq 
implies a planetary mass of msinz = 14.0 ± 0.8 Mj, placing the companion at 
the planet/brown dwarf boundary. Given its physical characteristics, the HD 
13189 system is potentially unique among planetary systems, and its continued 
investigation should provide invaluable data to extrasolar planetary research. 

Subject headings: stars:individual (HD 13189) — planetary systems — stars:early- 
type — stars:abundances — stars: atmospheres — stars: fundamental parameters 
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1. INTRODUCTION 

Ascertaining the physical properties of planetary host stars, the chemical abundances 
thereof in particular, is a critical component of understanding the formation and evolution 
of extrasolar planetary systems. Of the approximately 131 stars presently known to harbor 
planets, Fe abundances have been derived for no fewer than 117, with the abundances of 
numerous other elements having been determined for many. The results of these abundance 
analyses have led to the now well-known discovery that stars with planets tend to be metal- 
rich compared to random field stars (e.g., Fischer & Valenti 2005; Santos et al. 2005). Two 
hypotheses have emerged as possible explanations of the planet-metallicity correlation: ac- 
cretion of H-depleted, rocky material by planetary hosts and preferential planet formation in 
high-metallicity, proto-planetary disks. Both of these propositions are discussed extensively 
in the pertinent literature (e.g., Gonzalez 1997; Ida & Lin 2004; Fischer & Valenti 2005). 

An interesting consequence of the radial velocity (RV) method used currently to detect 
extrasolar planets is the limited range of stellar spectral types of stars chosen as targets. 
Planet searches generally focus on older, main sequence (MS) late-F, G, and K dwarfs, be- 
cause it is these stars that are bright enough to obtain high signal-to- noise, high-resolution 
spectra, have an ample number of usable spectral lines for the RV analyses, and have rota- 
tion rates and activity levels that facihtate the detection of planets. Furthermore, there is 
potentially great sociological significance in finding planets around stars like our Sun. While 
unique challenges exist for finding planets around younger stars and stars of earlier and later 
spectral types (e.g., Setiawan et al. 2003; Paulson, Cochran, & Hatzes 2004), the paucity of 
known planets in orbits around these "other" stars prohibits a full elucidation of extrasolar 
planet formation and evolution. 

Recognizing the importance of filling the gaps of the currently known planetary host 
sample, a handful of groups have initiated dedicated surveys of low mass stars (Endl et 
al. 2003) and evolved G- and K-giants, the progenitors of which were massive early-type 
MS dwarfs (Hatzes et al. 2003). Most recently, the Tautenburg Observatory Planet Search 
(TOPS) program, using the Alfred- Jensch 2-m telescope at the Thiiringer Landessternwarte 
Tautenburg (TLS) with follow-up observations with the 2.7-m Harlan J. Smith and Hobby 
Eberly (HET) telescopes at The McDonald Observatory, has announced the discovery of a 
giant planet orbiting the K2 II giant HD 13189 (Hatzes et al. 2005; Paper I). Based on an 
estimated luminosity of 3.6 Lq and the evolutionary tracks of Girardi et al. (1996), Paper 
I estimated a mass range of 2-7 Mq for the MS progenitor of HD 13189. Thus, HD 13189 
may be the most massive star known to harbor a planet. Here, we present the results of our 
abundance analysis of this potentially unique planetary host. 
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2. OBSERVATIONS AND DATA ANALYSIS 

Our abundance analysis makes use of two template spectra (spectra taken without an 
iodine cell) obtained as part of the TOPS program on 2002 December 19 using the Alfred- 
Jensch 2-m telescope and the high-resolution coude echelle spectrometer at TLS, located in 
Tautenburg, Germany. The spectrometer consists of an echelle grating with 31.6 g mm~^ 
and a cross-dispersing grism; the VISUAL grism has been used, along with a single 2048 x 
2048 CCD detector with 15 /jlir pixels, providing a wavelength coverage of 4660-7410 A. Each 
spectrum has a resolution of R — 67, 000 and a typical signal-to- noise (S/N) ratio of ~ 115. 
A standard data reduction process, which includes bias subtraction, flatfielding, scattered 
light removal, extraction, and wavelength calibration, using the usual routines within the 
I RAF facility has been applied to the spectra. 

The two individual HD 13189 spectra were co-added into a single, higher S/N ratio 
(S/N ^ 160) spectrum, and equivalent widths (EWs) were measured by fitting lines with 
Gaussian or Voigt profiles using the one-dimensional spectrum analysis package SPECTRE 
(Fitzpatrick & Sneden 1987). Spectral line selection was aided by the collection of Thevenin 
(1990); only lines labeled clS C3jS6 Oil '- lines with estimated internal uncertainties of < 0.05 dex 
in oscillator strengths (log (7/) and presumably free from blends in the solar spectrum- were 
considered. The hst was culled further by examining each fine in the spectrum of HD 13189 
and eliminating those that were deemed to be blended or otherwise unmeasurable. Our final 
line list, along with the measured EWs for both HD 13189 and the Sun, is presented in 
Table 1. Solar EWs have been measured in a high-quality {R = 60,000 and S/N = 1,000) 
spectrum of the daytime sky obtained with the 2.7-m Harlan J. Smith telescope at the 
McDonald Observatory (Schuler et al. 2005b). 

Stellar parameters have been determined spectroscopically via an iterative process. The 
Teff and ^ are determined by adjusting their values until there are no correlations between 
line- by- line [Fe/H]^ abundances (as derived from Fe I lines) and excitation potential (x) and 
reduced EW (logEW/A), respectively. This method for determining T^s and ^ is used often, 
but care must be taken to ensure that there is no ab initio correlation between excitation 
potential and reduced EW, which can lead to degenerate solutions. This is the case with 
our Fe I list, which consists of 86 fines. The list was trimmed until the correlation between 
excitation potential and reduced EW was efiminated; this left 47 Fe I fines with a satisfactory 
range of excitation potentials (0.91 < x ^ 4.99 eV). Uncertainties in Teff and ^ (la) are 
determined by adjusting the parameters until the corresponding correlation coefficient has a 
la significance. Finally, log g is fixed by forcing Fe abundances derived from Fe I and Fe II 
lines into agreement; the uncertainty in \ogg is based on the uncertainty in the difference of 
the Fe I and Fe II abundances. The final stellar parameters and uncertainties are given in 
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Table 2. We note our derived log 51 value supports the finding of Paper I that HD 13189 is 
a highly evolved giant. 

Abundances have been derived using the LTE stellar hne analysis package MOOG (Sne- 
den 1973; Sneden 2004, private communication). The abfind driver was used to force-fit 
abundances to the line-by-hne EWs for all elements except O. We have used the blends 
driver to derive the abundance of O from the A6300 [O I] line, taking care to account for 
the non- negligible contribution to the feature by a Ni I blend (e.g., Allende Prieto et al. 
2001). The specifics of our O analysis follows that of Schuler et al. (2005a), which should be 
consulted for details. Model atmospheres with the convective overshoot approximation for 
HD 13189 and the Sun have been interpolated from the ATLAS9 grids of Kurucz; the solar 
parameters — 5777 K, logg^ = 4.44, and ^ = 1.25 were adopted. Atomic parameters for 
all transitions are from Thevenin (1990), except for those of Mg, which are from the VALD 
database (Piskunov et al. 1995; Kupka et al. 1999; Ryabchikova et al. 1999), and that of 
A6300 [O I], which is from Allende Prieto et al. (2001). 

Final absolute abundances for HD 13189 and the Sun are given in Table 1, and the 
final stellar parameters, relative abundances, and uncertainties for HD 13189 are presented 
in Table 2. The Fe I abundances are derived with the complete fist of 86 fines. Relative 
abundances are the mean of a line-by-line comparison with solar values, thus limiting the 
impact of possibly inaccurate ^'/-values. Final abundance uncertainties are the quadratic 
sum of the uncertainties related to the adopted stellar parameters- based on sensitivities to 
arbitrary changes in T^s, logg, and ^- (Table 3)- and uncertainties in the mean abundances. 
For Mn and O abundances, which are based on a single feature, the final uncertainties 
incorporate uncertainties in the solar abundances. Sensitivities to the adopted C (important 
for molecular equilibrium calculations, primarily related to CO) and Ni abundances are also 
included in the final O abundance uncertainty. The final abundance uncertainties are on 
the order of 0.05-0.10 dex, but larger uncertainties due to systematic errors, e.g., related to 
relative differences in the atmospheric models, cannot be ruled out. 



3. DISCUSSION 

The derived Fe abundance, [Fe/H] = — 0.58±0.04, places HD 13189 in the lower envelope 
of the planetary host metaUicity distribution and makes HD 13189 one of the most metal- 
poor planetary hosts yet discovered. Furthermore, with a semi-major axis in the range 



^The bracket notation is used to denote abundances relative to solar values, e.g., [Fe/H] = 
log{N{Fe)/N{H)}^ - \og{N{Fe)/N{H)}Q, where logiV(if) = 12.0. 
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1.5-2.2 AU (Paper I), the planetary companion of HD 13189 is the only one known to have 
a semi-major axis > 1.0 AU and orbit a star with a metallicity < —0.40, a combination 
that is potentially important to understanding giant planet formation (Rice & Armitage 
2005). Abundances of other elements considered here are sub-solar, and the [m/Fe] ratios are 
unspectacular when compared to the general field population (e.g., Edvardsson et al. 1993; 
Fulbright 2002). The planet-metallicity correlation is one of the most scrutinized results to 
emerge from extrasolar planet studies. Fischer & Valenti (2005) offer the most thorough 
investigation of the correlation to date; their conclusion, based on a uniform analysis of 
850 stars with and without planets, is that gas giant planets preferentially form in high- 
met allicity environments. Various groups have also argued for this "primordial" scenario 
(e.g., Pinsonneault, DePoy, & Coffee 2001; Santos 2003), and there are two main pieces of 
evidence given in support. First, there is no correlation between convection zone (CZ) depth 
and metallicity of planetary hosts, an effect that is expected if accretion of H-depleted rocky 
material has occurred. Second, the metallicity distribution of planet-harboring subgiant 
stars does not differ from that of dwarfs; the deepening CZ of a subgiant should dilute any 
accretion signatures that might have occurred. 

Despite the substantial evidence pointing to the primordial scenario for planetary for- 
mation, data in support of the accretion of H-depleted rocky material hypothesis also exist, 
albeit they are more sparse. Gonzalez (1997) was the first to report on the super-solar 
metaUicities of planetary hosts and to suggest the correlation may be a result of accretion. 
Shortly thereafter, Israelian et al. (2001) detected ^Li in the metal-rich planetary host HD 
82943; ^Li is a volatile species that should be destroyed beyond the level of observabihty in 
the atmospheres of solar- type stars during pre-MS evolution (e.g., Proffitt & Michaud 1989). 
Israelian et al. interpreted the presence of ®Li- which should be preserved in giant planets- 
in the HD 82943 as possible evidence of accretion; the possibility of heightened ^Li abun- 
dances in the atmospheres of stars that have engulfed rocky material was later confirmed by 
Montalban & Rebolo (2002). The ^Li detection was subsequently challenged by Reddy et 
al. (2002), who questioned the atomic data used by Israelian et al. In response, (Israehan 
et al. 2003) reanalyzed ^Li in HD 82943 using updated atomic data and confirmed their 
original detection. Further evidence for the accretion hypothesis comes from the study of 
Smith, Cunha, & Lazzaro (2001), who found a trend of increasing abundance for elements 
with increasing condensation temperatures (Tc) for a small subsample of planetary hosts. If 
accretion of rocky material has occurred, one might expect the material to be rich in high-Tc, 
refractory elements compared to low-Tc, volatile elements (Gonzalez 1997). However, simi- 
lar Tc-dependent abundance trends in planetary hosts have not been found by others (e.g., 
Sadakane et al. 2002). 

Additional support for the accretion hypothesis may be emerging from the discovery and 
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analyses of giant stars with planetary companions. As of this writing, planets are known 
to orbit two G-type and five K-type giants (including HD 13189), and their metallicity 
distribution is not similar to that of planet-hosting dwarfs. In Figure 1 we plot relative Fe 
abundances versus Tcfj for the seven giants (red circles) and for planet hosting dwarfs from 
Fischer & Valenti (2005; black circles), from Santos, Israelian, & Mayor (2004; black squares), 
excluding duplicates with Fischer & Valenti, and from Santos et al. (2005; black triangles), 
again excluding duphcates. The horizontal hne represents the mean [Fe/H] abundance of 
the combined, 126 dwarf sample. The giant data for HD 47536, HD 59686, HD 137759, and 
HD 219449 are from Sadakane et al. (2005), and data for HD 11977 and HD 104985 are from 
Setiawan et al. (2005) and Sato et al. (2003), respectively. In the combined dwarf sample, 
41.3% have derived [Fc/H] abundances higher than +0.20 dex, and if the giants follow the 
same metallicity distribution, we would expect ~ 3 to have [Fe/H] abundances higher than 
+0.20 dex (see also Sadakane et al. 2005), which is clearly not seen. If the dwarf and 
giant samples are drawn from the same parent population, there is only a 0.6% probability 
of obtaining the observed metallicity distributions according to a two-sided Kolmogorov- 
Smirnov test. 

While the sample of giant stars with planetary companions is not statistically signifi- 
cant, the indication of a metallicity distribution that possibly differs from that of dwarfs is 
intriguing. One possible explanation for the difference is that the planet-metallicity corre- 
lation found for the dwarfs is a result of accretion, and the giants have sufficiently diluted 
their atmospheres via evolution-induced mixing. This conclusion, however, would be difficult 
to reconcile with current results for subgiants (e.g., Fischer & Valenti 2005). Alternatively, 
planet formation may be sensitive to the mass of the host star, independent of metallicity. 
The progenitors of G and K giants are early-type MS dwarfs and would be more massive 
than the F, G, and K dwarfs known to have planets. Thus, the current results for giants 
with planets would support a mass- dependent planetary formation scenario. Furthermore, 
Fischer & Valenti showed that the occurrence of planets around stars increases linearly with 
stellar mass. However, they also found a correlation of rising metallicity with increasing 
planetary host mass and argued that the increased occurrence of planets around more mas- 
sive stars is most likely spurious. A larger number planets orbiting subgiant and giant stars 
need to be analyzed in order to determine if either of these explanations are plausible, or 
if the difference in metallicity distributions is due to other effects, such as different plane- 
tary formation mechanisms (e.g., Chauvin et al. 2005). Including more evolved stars in RV 
surveys is highly encouraged. 

Our analysis provides the opportunity to place further constraints on the mass of HD 
13189 from Paper I. The tracks of Girardi et al. (2000) are plotted in the Teg-logg plane in 
Figure 2; the position of HD 13189 is marked by the open star. From our derived stellar 
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parameters, we find a mass of M^, — 3.5 M© for HD 13189. Unfortunately, the range of 
possible masses, = 2-6 M©, is not well-constrained due to the uncertainty in the logs' 
value and is not a significant improvement over that given by Paper I. Adopting a stellar 
mass of 3.5 M©, we find a minimum companion mass of 14.0±0.8 Mj, placing the companion 
at the planet/brown dwarf boundary. HD 13189 and its companion are a potentially unique 
system among planetary systems known today. The star has one of the lowest metallicities 
and possibly the highest mass of all planetary hosts, and the companion likely lies on the 
planet/brown dwarf boundary and has the largest semi-major axis for a planet orbiting a 
significantly metal-poor star. Further investigation of the HD 13189 system is needed to place 
more stringent constraints on the mass of the two objects and should provide unparalleled 
data in the effort to understand the formation and evolution of planetary systems. 

REU students J.H.K. and M.T. gratefully acknowledge the support of summer research 
through the NSF-REU Award 0353849 to Clemson University S.C.S. and J.R.K acknowledge 
support for this work by grant AST 02-39518 to J.R.K. from NSF, as well as a generous grant 
from the Charles Curry Foundation to Clemson University. 
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Table 1. Equivalent Widths AND Abundances 





A 


X 




EWo 




EW 




Species 


(A) 


(eV) 




(mA) 


logTVe 


(mA) 


log AT 


Fe I . . . 


5633.95 


4.99 


-0.27 


82.4 


7.59 


95.1 


6.97 




5705.47 


4.30 


-1.57 


40.1 


7.60 


82.7 


7.14 




5720.90 


4.55 


-1.95 


20.1 


7.77 


46.4 


7.30 




5731.77 


4.26 


-1.19 


59.7 


7.54 


97.1 


6.92 




5732.30 


4.99 


-1.50 


16.9 


7.63 


31.7 


7.16 



Note. — Table 1 is presented in its entirety in the electronic edition of the 
Astrophysical Journal. A portion is provided here for guidance regarding its 
form and content. 
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Table 2. Final Results 
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Teff (K) . . 
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0.22 
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-0.58 


0.04 


[0/H] .... 
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0.08 


[Na/H] . . . 


-0.31 


0.06 


[Mg/H] .. 
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0.05 
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0.06 


[Si/H] .... 
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0.07 


[Ca/H] . . . 
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0.07 
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0.07 
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-0.42 


0.07 


[V/H] .... 


-0.38 


0.12 


[Cr/H] . . . 


-0.47 


0.09 


[Mn/H] . . 


-0.58 


0.09 


[Ni/H] . . . 


-0.57 


0.04 
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Table 3. Abundance Sensitivities 
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01 ... 
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Al I . . . 


±0.13 


±0.00 
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Si I . . . 
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Cal .. 


±0.18 


0.00 


^0.10 


Sc II . . 


:f0.04 


±0.11 


^0.06 


Ti I . . . 


±0.23 


0.00 


:F0.09 


VI.... 


±0.26 


±0.02 


^0.15 


Cr I... 


±0.12 


±0.01 


^0.08 


Mn I. . 


±0.09 


±0.03 


TO.02 


Ni I . . . 


±0.02 


±0.07 


^0.10 
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Fig. 1. — Relative Fe abundances of giants (red circles) and dwarfs with planetary compan- 
ions versus T^s. The sources of the giant abundances and Teg other than those for HD 13189 
are given in the text. Dwarf abundances are from Fischer & Valenti (2005; black circles), 
Santos, Israelian, & Mayor (2004; black squares), and Santos et al. (2005; black triangles). 
The horizontal line represents the mean Fe abundance ([Fe/H] = +0.13) of the dwarf sample. 
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Fig. 2. — Evolutionary tracks of Girardi et al. (2000) for (from top left to bottom right) = 
1.5, 2.0, 3.0, 3.5, 4.0, 5.0, 6.0, and 7.0 Mq. The tracks are characterized by a metalhcity 
of [Fe/H] = —0.50 or Z ^ 0.008. The open star with error bars marks the location of HD 
13189 in the Teg-log (7 plane. 



